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KATHRIN BRINGMANN 1 AND JAN MANSCHOT 

Abstract. We determine asymptotic formulas for the coefficients of a natural class of 
negative index and negative weight Jacobi forms. These coefficients can be viewed as a 
refinement of the numbers pk(n) of partitions of n into k colors. Part of the motivation for 
this work is that they are equal to the Betti numbers of the Hilbert scheme of points on 
an algebraic surface S and appear also as counts of Bogomolny-Prasad-Sommerfield (BPS) 
states in physics. 
< 

<N 

[_h \ 1. Introduction and Statement of Results 

Jacobi forms were first systematically studied by Eichler and Zagier [T3] and enjoy a wide 
variety of applications in the theory of modular forms, combinatorics [231 EH] , conformal field 
theory [IHI2T], black hole physics [HIES], Hilbert schemes of points JTH], Donaldson invariants 
[PT] . and many other topics. This paper focuses on a class of negative index Jacobi forms 
with a single order pole in the elliptic variable w. The analysis of the coefficients of such 
functions is more complicated then the well-understood class of Jacobi forms which depend 
holomorphically on w. It turns out that these Fourier coefficients (in w) are not modular 
but related to quantum modular forms [HI [36] . The appearance of these functions in the 



above mentioned topics calls for an explicit knowledge of their coefficients and in particular 
of their asymptotic growth. In this paper we provide such asymptotic formulas. One of the 
immediate motivations is the counting of BPS states in physics and in particular those with 
vanishing angular momentum. This motivation is explained in more detail after stating the 
results. 

1.1. Statement of results. We consider the following class of negative weight 1 — k/2 
d', (k £ N, k > 3) and index —1/2 Jacobi forms 

h k (w;r) :-- 



8i(w;r)ri(T) k 3 ' 
with (q := exp(27rir), ( := exp(2niw)) 

n>l 



n>l 
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We are interested in two expansions of hk{w',r). The first expansion is in terms of the 
coefficients a m ,k(n), defined by: 

(1-1) 9* (C* " C h ) h k (w; r) =: ]T a m , k (n)( m q n , |C<?I, \C\\ < 1- 



n>0 

mgZ 



Note that ^ m£Z a rai )t(n) = Pk(n), where pk(n) denotes the number of partitions of n into k 
colors. These are enumerated by rj(r)^ k : 

^2Pk(n)q n ~^ 



77 (Y 
n>0 /v 



,fc' 



The second expansion is motivated from physics, and is based on the fact that the coeffi- 
cients of q are Laurent polynomials, symmetric under ( — > £ -1 and with maximal degree n. 
Therefore, we can express h k (w; r) as: 

(1-2) q^ (d - C~*) h k (w; r) = J^ 6 m , fc (n) X2m+i (C*) Q n , 

m,n>0 

with 

Xm(C) : = 737~i ' b m: k(n) := a m ,k{n) - a m+ x jk (n). 



Following the approach of Wright [34J, we determine the asymptotics of the coefficients 
o>m,k(n) for large n. 

Theorem 1.1. We have for N > 1: 

N / n ox-\ 1+ ^+l 

/ \ /^ \- h S~^ r I n\ V^ / -\ 2+2l+fc I TTy/K — Z4j \ 

a m ,k{n) = (2tt) 2 2^«m,fc(«) 2^ PfcU) n 4 I 7^ J 

(I Ji \ I iV+3 k TrV3kn \ 

tc\ j~{k-2Aj)n\ +0 In 2 i e 2v ^ I, 

where d m> k{i) are defined by equation ( PQ| ) ; Pk{j) denotes the number of partitions of j into 
k colors, and Ii(x) is the usual I -Bessel function. 

Theorem 11.11 allows us to compare the asymptotic growths of a m ^ k (n) for different values of 
m. The asymptotic behavior of the Bessel function : 



e x 



h{x) = -== (1 + O (x- 1 )) 



directly yields: 
Corollary 1.1. We have 



fc+7 
■' :: I k\ 4 /2fc7~ / _ l-l I ,'c /2fc77 



.. .. "3/9 9\/ ^Jj-K / flj \ / Zkn I Id-He / ^fc'. 

Om, fc (n) - o r ,ife(n) = ?r (r 2 - m 2 ) (8n) 4 _ e^v 3 + O [n 4 e ^v 3 

From Theorem 11.11 the asymptotics of b m ^ k (n) for large n immediately follow: 
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fc+7 



Corollary 1.2. We have 

(1.3) 6 m , fc (n) = (2m + l)7r 3 (8n)"— (-J e ff V— + (n-^e'V— J 

Note that this corollary shows that b m k(n) increases with m in the limit of large n. Beyond 
the validity of equation (II. 3p . b m ^{n) eventually decreases with increasing m for fixed n, and 
in particular b m ^{n) = for m > n. 

We next compare the asymptotic behavior of the coefficients bo,k{n) with those of Pk{n). 
It is well-known that the asymptotic growth of the latter is given by [TBI [27] : 



k\ 4 



3+fc „ / 2kn ^ ( 5 + fc „ I 2kn 
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p fc (n) = 2[-l (8n)- — e 7T v—+0(n- — e T: 
Thus we find for the ratio bo } k{n)/pk{n): 

bo, k (n) 7T 3 ( kY 



Pk{n) 



16 \3nJ V ; 



1.2. Motivation: Moduli spaces and BPS states. BPS states of both gauge theory 
and gravity have been extensively studied in the past for a variety of reasons. These states 
are representations of the SU(2) spin massive little group in four dimensions labeled by their 
angular momentum or highest weight J. The subset of BPS states with vanishing angular 
momentum (J = 0), also known as "pure Higgs states" [4], have recently attracted much 
interest [4J [25J ESI ED]- The states with J = are in some sense more fundamental. In 
particular in gravity, these states are candidates for microstates of single center black holes, 
and as such are the states relevant for studies of the Bekenstein-Hawking area law of black 
hole entropy. 

String theory relates BPS states to the cohomology of moduli spaces of sheaves supported 
on a Calabi-Yau manifold. From this perspective the SU(2) spin representations correspond 
to representations of the Lefshetz SL(2) action on the cohomology of the moduli space [32J. 
The states with J = correspond to the part of the middle cohomology which is invariant 
under the Lefshetz action. 

In the present work, we consider moduli spaces of semi-stable sheaves supported on a 
complex algebraic surface S, which can be thought of as being embedded inside a Calabi- 
Yau manifold. If S is one of the rational surfaces, the sheaves can be related to monopole 
or monopole strings in respectively four and five dimensional supersymmetric gauge theory 
through geometric engineering [20| [26] . If S is a K3 surface, the sheaves correspond to 
(small) black holes in M = 4 supergravity also known as Dabholkar-Harvey states [7J [8] . We 
specialize to the moduli space of sheaves with rank r = 1 and 1st and 2nd Chern classes 
c\ G H 2 (S, Z) and c 2 G H A (S, Z). These moduli spaces are isomorphic to the Hilbert scheme 
of C2 points on S (viewing ci as a number). Gottsche has determined the generating function 
of the Betti numbers of the Hilbert schemes [16J. We need to introduce some notation to 
explain his result. 
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Let A4(n) be the Hilbert scheme of n points. Let furthermore P(X, () := X}i=o @i{X) C 

be the Poincare polynomial of X with fii(X) the ith Betti number of X. We choose the sur- 
face S such that (3i(S) = {3s(S) = 0. Then, we have [16J: 



n>0 



■^ M( rip(M(n),C)q n = q Sai *F 1 fa* - C~* 



e^w-T)^)^)- 1 ' 



which precisely equals the function in ( 11. ip with k = ^(S) + 2. The coefficients a m ^{n) 
are in this context the Betti numbers of the moduli spaces. The expansion (II. 2p in terms 
of bj t k(n) decomposes the cohomology in terms of (2J + 1) -dimensional SL(2) or SU(2) spin 
representations. For J = and k = 24, our formula (11.31) confirms nicely the numerical 
estimates for 6o,24(^) m [HI Appendix C]. Note that since this analysis is carried out in the 
so-called weak coupling or D-brane regime, and the coefficients bj t ^(n) are not BPS indices, 
these coefficients can not be claimed to count black holes with fixed angular momentum. 

In this example of the Hilbert schemes of points, we find a few qualitative differences 
between the asymptotics of the number of pure Higgs states and the remaining BPS states 
compared to [UHO]. In 0], it is shown that the number of pure Higgs, bo,k(n), is exponentially 
larger than the number of remaining BPS states given by p m ,k{n) — bo,k(n). On the other 
hand, equation (11.31) shows that in the context of this paper, the number of pure Higgs 
states and the other states have the same exponential growth but differ in the power of n. 
Moreover, the number of SU(2) spin multiplets increases with J for small J. Of course, we 
could not expect to find similar qualitative features since the limit for the asymptotics are 
of a different nature. 

To make this more explicit, note that in [U [10] moduli spaces of quiver representations 
are considered, and that sheaves on toric surfaces also allow a description in terms of quivers 
(now with a potential) [6j [29]. Whereas in [U [10] a limit is taken with a large number of 
arrows but the dimension of the spaces associated to the nodes kept fixed, the Hilbert scheme 
of n ^> 1 points corresponds to increasing dimensions of the spaces associated the nodes, 
with the number of arrows kept fixed. It will be interesting to extend the techniques in the 
present paper to partition functions for higher rank sheaves on surfaces [231 EU HH], and 
partition functions of black holes and quantum geometry [5J [H)J [IS] • 

The paper is organized as follows: In Section [2] we rewrite the functions of interest in 
terms of false theta functions and determine their Taylor expansion. Section [3] uses the 
Circle Method to prove our main theorem. 
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2. Relation to false theta functions 

We start by writing the generating function of a mt k{n) for fixed m and k in terms of the 
functions i9 m (g) defined by: 

d m (q) := (1 + q m ) ^(-1)^^+™ - 1. 

n>0 
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These functions are examples of false theta functions, which were first introduced by Rogers 
[28] and have attracted a lot of interest recently. Using the Rogers and Fine identity one 
can relate d m to so-called quantum modular forms which are functions mimicking modular 
behavior on (subsets of) (Q 



Proposition 2.1. We have for |£g|, \( 1 q\ < 1: 

withiqU^nZiii-q")- 

Proof: We first note that the inverse theta function 9\(w\ r) _1 can be written as a sum over 
its poles [3]: 



iql ^2 — C 2 J 1 



{-l) n q 



n(n + l) 



M -^^-WJ-OT^ 



6 1 (w;r) (g)s/ *' *£ 1 - (q* 

Using geometric series expansion, we may rewrite (12.11) as 

» V — "V , < n(n + l) . / -1 \ V — "V , ^ n(n-\-l) . 

ttf(i - o j^i- 1 ) <i* + c + 7y(!- c 1 ) j^i- 1 ) q> c m - 

[Q)oo n > WJoo n > 

m>0 m>0 

From this the statement of the proposition easily follows. □ 

The function d m is not modular but may be nicely approximated by its Taylor expansion. 
For this we use the following general lemma (see [35] for the case of real functions). 

Lemma 2.2. Let f : C — > C be a C°° function. Furthermore, we require that f(x) and all 
its derivatives are of rapid decay for Re(x) — > oo. Then for t —$■ with Re(t) > and a > 0, 
we have for any iV 6 N : 

£ / «m + a)f) = \r md X - ± ^^e + O («»«) . 

m=0 Jo n=0 ' 

Here B„(x) denotes the nth Bernoulli polynomial. 

To use Lemma 12.21 we write for fixed N > 1 and q = e~ z 

N 

(2.2) q^ m {q) =: £ <W(V + O {z N+1 ) . 

e.=i 

Lemma 12.21 then gives 

Lemma 2.3. We have for N > 1 

N 
$m (q) = (1 + q m ) g~K™+l) 2 J2 Cm{ey -1 + (^ +1 ) 

1=1 

with 

, n . {-ly-^ ( fm 1\ (m 3 

Cm{i) '= *!(2*+l) ^ 2£+1 U + 4 ) ~ 2l+1 U + 4 



In particular the first values for d mj k(£) are: 

mA) = 4' mA) = ~96 + 16' mA) = ~ 16 " 384 + 4608 + 192' 

Proof: We may write ^m(Q) as: 

^) = (i +9 "-)^(™ + ^E(/(("+f + i)^)-/(("+f4)^))-i 

with fix) := e~ 2a: . Substitution of Lemma E21 gives the desired result. □ 

Remark. We note that the case m = can be easily concluded from [34] . 



3. Use of the Circle Method 

In this section, we prove Theorem II .11 following an approach by Wright [51j . To prove the 
theorem, we assume via symmetry that m > and set 



F m ,k{q) ■= ^Glm,fcW<f 



n>0 



By Proposition 12.11 we obtain that 



F m ,k(q) = p^m(?)- 



'OG 



By Cauchy's Theorem, we have for n > 1 

where C is a circle surrounding counterclockwise. We choose e~ v for the radius of C with 
77 = 7tyjr- and split C into two arcs C = Ci + C2, where Ci is the arc going counterclockwise 
from phase —2-q to 2rj and C2 is its complement in C. Consequently, we have 

a mtk {n)=M + E 

with 

£ := — / ^# rfg . 

We will show that the main asymptotic contribution comes from M. Moreover we parame- 
trize q = e~ z with Ke(z) = 77. 
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3.1. The integral along C±. In the integral along &_, we approximate F m ^ by simpler 
functions. Firstly, recall that from the transformation law of the 7/-function we obtain [2] 

, s l Fz z ,^ i 

(3.1) — = \— e"53 + ^ 

(e-^e-*)^ V27T i -'- -'- 



Thus we want to approximate ^-L- by 

W/oo 

^2(27r) _ 2e _ 2fP fc (e~ 



i+ x) p*(^U" 



" 24 



\ 24l-k<0 



where 

/ 

To be more precise, we split 

M = Mi + Ex 
with 

We first bound i?i which turns into the error term. Firstly we obtain from (13.11) 

i h. 



— e-wR, e -— =0(1). 



(e-*;e-*)« V27r 
To bound $ m , we use that on Ci 

|^| 2 = ?7 2 + Im(^) 2 < rf + 4?7 2 . 



Thus, by Lemma [2.31 

\& m {q)\ < M <?7, 
where throughout g>(x) <C f(x) has the same meaning as g(x) = 0(f(x)). Using that the 
length of C\ is 0(7/), we may thus bound 

-Ei < n ^ l e ^ . 

We next investigate M\. We aim to approximate d m by its Taylor expansion given in 
Lemma 12.31 and thus we split 

Mi = M 2 + E 2 

with 

i N r i * 2 



/ Cl g™+! V2vr, 



h - / 



e=i 



We first estimate E 2 and show that it contributes to the error term. By Lemma 12.31 



|£ 2 |< / e n "|zp 

>Cl 



- N ^e^ R < l *)dz. 



Since 



we may bound 



Moreover on C% 



Re(z) 1 



\z\ 2 Re(z) r] 



ir 2 k„ /1\ 2kn 

nr]+ — Re U )- n y-r 



x 2 + y 2 < rj 2 + 4?7 2 < T] 2 . 



As before, the path of integration may be estimated against 77. Thus 

\E 2 \ < r ? Ar+ 2+ 2 e 7r V— < n-T-i-ieV — . 



We next decompose 

JV 

(3-2) M 2 = (2vr)-f^rf m , fc (£) £ ^(j)^., 

£=1 24j-fc<0 

where for s > 



T (fe) . 



2VT2 ./ c 



z e $z 



'- + (n+l)z 



dq. 



These integrals may now be written in terms of the classical J-Bessel functions. 
Lemma 3.1. We have 



X sJ = n 2 



s+l 



lfn Vk ^ 2Aj y^I-s-iln ] ll(k-24j)n)+0[n^-ie^r 



Proof: Let T> be the rectangular counterclockwise path from —00 — 2it] to — 00 + 2ii] with 
endpoints r\ — 2ir\ and r\ + 2ir\. Denote by V iy i = 1,2, 3, the paths 1) from —00 — 2ir\ to 
77 — 2ir], 2) from 77 — 2irj to 77 + 2z?7, and 3) from 77 + 2^77 to —00 + 2%r\. Making the change of 
variables q = e~ z gives that 

1 



X 



(k) 



2-ni 



z e sz 



*'dz. 



v< 



We next use the Residue Theorem to turn this integral into an integral over T>. For this we 
bound the integrals along T>\ and V3. We only give the details for V3. On this path we may 
bound 



Re 



Writing z = 77(1 + 2i) — u, < u < 00 gives 



1 1 

< — < — • 

\z\ 2rj 



(77 — u) + Arf <C 77 + u. 
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Thus the integral along V 3 may be bounded by 

(3.3) < e 3 ^(» _i ) / (77 + u) s e n ^- u) du < e^T +n " ( 77 s / ' e- nu dw + / w s e- nu tk ) . 



'v 

The second term is an incomplete Gamma function and thus exponentially small. Thus (13. 3p 
may up to an exponentially small error be bounded by 



TV k 1 \ , - -1 1 7ry/3kn 

rf^+^L (1 _ e -m} < n -i-f e ^73-. 



s 7r\/3fc-n 

(l-e~"'') <n _i ~ 

n 

In the remaining integral we make the change of variables z — - to get 



(3.4) Xf/ = n - * -1 — / tV + — — dt + O n^-ie ^ 

2-7U J j, \ 

We now use the following representation of the J-Bessel function pQ 

'(0+) 



I e (2y/z) =z*^-. [ r^exp (t + ~) dt, 
2m J_ 00 \ t) 



where the integral is along any path looping from —00 around back to —00 counterclock- 
wise. Substitution into ( 13. 4 p gives the claim. □ 

Substitution of Lemma 13.11 in equation ( 13.21) yields 

-\\p 1 h. / 1 ' 

M 2 = ( 27 r)-f J2 ^ffl^(j>' g±fM ( 7rVfc 7 g 24j ) '/_«_$ (*^(*-24j)n 

24j-fc<0 

/ _3_fc 7iV3fcn 
+ O I U 2 4 e 2^2 

3.2. The integral along C 2 - On C 2 , Im(z) varies from — 2i] to — 2n + 2i]. Using a rough 
bound for the theta function, we find 

|tf m (g)| < 2V e -f(-+l+2 m )Re(,) + 1 <<2 ^ e -RoW = _J_ < i 
n>l n>0 ' 



Using ( 13. ip we obtain the bound 



1 



Refi' 



(3.5) — < e 

Now 

/l\ _ 7] 1 

\2/ 77 2 + \m.(z) 2 ~ 5?7 
Thus ( 13. 5p may be estimated against exp(^-). This gives that: 

E < eW+"" < e f ^ 

This is exponentially smaller than the other errors. Combining this with the results of 
Subsection 13.11 therefore gives Theorem 11.11 
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